Weak Turbulent Kolmogorov Spectrum for Surface Gravity Waves 
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We study the long-time evolution of surface gravity waves on deep water excited by a stochastic 
external force concentrated in moderately small wave numbers. We numerically implemented the 
primitive Euler equations for the potential flow of an ideal fluid with free surface written in Hamil- 
tonian canonical variables, using the expansion of the Hamiltonian in powers of nonlinearity of up 
to terms of fourth order. 

We show that due to nonlinear interaction processes a stationary Fourier-spectrum of a surface 
elevation close to < |77fcp >~ fc"^/^ is formed. The observed spectrum can be interpreted as a 
weak-turbulent Kolmogorov spectrum for a direct cascade of energy. 



Kolmogorov was born in 1903. Now, in the year of 
his centenary, his greatness is obvious not only for pure 
and applied mathematicians. All physicists appreciate 
his pioneering works on powerlike cascade spectra in tur- 
bulence of the incompressible fluid 0. It is obvious 
now that cascade processes, similar to the Kolmogorov 
cascade of energy, play a very important role in many 
different fields of physics, such as nonlinear optics [21, 
plasma physics 's^, hydrodynamics of superfluid He^, and 
so forth. 

In all these cases the physical situations are similar. 
There is an ensemble of slowly decaying, weakly nonlinear 
waves in a medium with dispersion. Such systems have to 
be described statistically. However this is not traditional 
statistical mechanics, because the ensembles are very far 
from thermodynamic equilibrium. Nevertheless, one can 
develop a systematic approach for the statistical study of 
weakly nonlinear waves. This is the theory of weak (or 
wave) turbulence . The main tools here are the kinetic 
equations for squared wave amplitudes. These equations 
describe the nonlinear resonant interaction processes tak- 
ing place in the wave systems. Like in the turbulence in 
incompressible fluid, these processes lead to the cascades 
of some constants of motion (energy, wave action, mo- 
mentum etc.) along the fc— space. In isotropic systems it 
might be either a direct cascade of energy from small to 
large wave numbers or an inverse cascade of wave action 
to small wave numbers 0|. In an anisotropic system the 
situation could be much more complicated |^. 

The brilliant conjecture of Kolmogorov still is a hy- 
pothesis, supported by ample experimental evidence. On 
the contrary, the existence of powerlike Kolmogorov spec- 
tra, describing cascades in weak turbulence, is a rigor- 
ous mathematical fact. These spectra are the exact so- 
lutions of the stationary homogeneous kinetic equation. 



completely different from the thermodynamic Rayleigh- 
Jeans solutions. 

Nevertheless, the case is not closed. The weak tur- 
bulent theory itself is based on some assumptions, like 
phase stochasticity and the absence of coherent struc- 
tures. This is the reason why justification of weak tur- 
bulent theory is an urgent and important problem. 

This justification can be done by a direct numerical so- 
lution of the primitive dynamic equation describing the 
wave ensemble. In pioneering works by Majda, McLaugh- 
lin and Tabak it was done for the 1-D wave system. 
The results obtained by these authors are not easily in- 
terpreted. In some cases they demonstrate Kolmogorov- 
type spectra, in other cases — power spectra with essen- 
tially different exponents. 

In article Q deviation from weak turbulent theory was 
explained by the role of coherent structures (solitons, 
quasi-solitons and collapses). If a 1-D system is free from 
coherent structures, weak-turbulent spectra are observed 
with a good deal of evidence |^, Hol ll]| . 

In spite of their heuristic value, the 1-D models so 
far developed have no direct physical application. Real 
physical systems, where wave turbulence is realized, are 
at least 2-dimensional. The most natural and impor- 
tant examples are capillary and gravity waves on deep 
water. A weak-turbulent theory of capillary waves was 
developed by Zakharov and Filonenko in 1967 JJj, who 
found that the correlation function of elevation r]{r, t) 
has to be < jr^^p >~ k''^^/^. This result was supported 
by laboratory experiment s, p erformed independently by 
three groups (in UCLA O, Niels Bohr Institute [l|| 
and the Solid State Physics Institute in Chernogolovka, 
Russia 15, 16J). The spectrum k~^^/'^ was obtained by 
a direct numerical simulation of Euler equation for in- 
compressible fluid with free surface by Pushkarev and 
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Zakharov [T1[T1[T3. 

The most interesting example of 2-D wave ensembles 
demonstrating weak-turbulent cascades is a system of 
gravity waves on the surface of deep water. We are sure 
that a weak-turbulent theory of these waves is key to un- 
derstanding the processes in a wind-driven sea. However, 
we do not concentrate on this point in our article. 

Our initial goal was to reproduce (and emulate), for 
gravity waves, the work which was done by Pushkarev 
and Zakharov 0] for capillary waves. One has to expect 
that this is a more difficult problem, because the lead- 
ing process in capillary waves is a three-wave interaction 
(dispersion of waves is of decay type), while for gravity 
waves the lowest order process is four-wave interaction 
(dispersion of waves is of non decay type). 

Attempts to perform direct numerical simulations of 
potential flow in an ideal fluid with a free surface 
were made by several authors |2^. Only in one arti- 
cle authors did pay interest to Kolmogorov-type weak- 
turbulent spectra [2l|. Authors of this paper observed 
the formation of Kolmogorov tails during the time evolu- 
tion of an artificially cut-off JONSWAP (Joint North Sea 
Wave Project, described in J2]) energy spectrum. The 
results of present article agree with the results of ji^l 
completely. However we would like to stress a difference. 

In our work, in contrast to , we study a forced tur- 
bulence, excited by external sources, posed in moderately 
low wave numbers. We show that the growth of wave en- 
ergy due to this forcing is arrested by nonlinear resonant 
four-wave processes, which leads to the formation of a 
powerlike Kolmogorov spectrum in the inertial interval. 
In this sense our article is a direct numerical confirmation 
of the weak-turbulent theory for surface gravity waves. 

Theoretical Background — We study the potential 
flow of an ideal inviscid incompressible fluid with the ve- 
locity potential ((> — (t>{x, y, z; t) 

in the infinitely deep domain occupied by the fluid. Equa- 
tions for the boundary conditions at the surface are the 
following 
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= 0. 



(1) 



Here ri{x, y; t) is the surface elevation with respect to still 
water, g is t he g ravity acceleration. Equations are 
Hamiltonian |22| with the canonical variables ri{x,y;t) 
and tp{x, y- 1) = 0(a;, y, t]{x, y; t);t) 



dri _5H dij _ SH 
dt S'tp ' dt Srj^ 
where H is the Hamiltonian of the system 
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dxdy grj 



Unfortunately H cannot be written in the close form as 
a functional of ry and ^p. However one can limit Hamilto- 
nian by first three terms of powers of rj and ijj 



H = H„ 



Hi 



Hi + H2 - 

^ ' ipktp ) dxdy, 
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dxdy. 



(3) 



H2 = - I ri{kijj) k{r/{kip)) + rjAip dxdy. 



Here A; is a linear integral operator (j: = \/~Aj , such 
that in fc-space it corresponds to multiplication of Fourier 
harmonics {ipj: = 57 / ijpe^^^dxdy) by ^k^ + fc^ . For 
gravity waves this reduced Hamiltonian describes four- 
wave interaction. Then dynamical equations Q acquire 
the form 



kip - {V{r]Vip)) - k[r]kip] 
+k{r]k[T]k'i{j]) + iA[7?2fci/,; 

(VV')2 - (fcV)2 
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-[kip]k[r]kil)] - [r]ki)]Aij + Dp+ Fp. 



(4) 



Here Dp is a dissipation term, which consists of hypervis- 
cosity on small scales and "artificial" damping on large 
scales. Fp is the driving term which simulates pump- 
ing on large scales (for example, due to wind). In the 
fc-space supports of and Fj: are separated by the in- 
ertial interval, where the Kolmogorov-type solution can 
be recognized. 

Another approach based on boundary integral approx- 
imation has been studied by Clamond and Grue in 24] . 

We study numerically the quasi stationary solution of 
equations 10}. According to the theory of weak turbu- 
lence, the Fourier spectrum of the surface elevation av- 
eraged by ensemble (corresponding to the fiux of energy 
from large scales to small scales) is 



< \Vkf >-- 



Cg^/^P^/^ 



(5) 



Here P is the energy flux, and C is a dimensionless Kol- 
mogorov constant. This spectrum is a stationary solution 
of the wave kinetic equation 12] . 

It is worth to say that an alternative spectrum was 
proposed earlier by Phillips |25j . In contrast to weak tur- 
bulent theory Phillips power-like spectrum can be formed 
due to a wave breaking mechanism and leads to a surface 
elevation spectrum < \rik\'^ >~ k~'^. It can be realized 
for higher level of turbulence, which is far beyond the 
weak turbulent theory. 

Numerical Simulation — For numerical integration of 
Q we used the following pumping and dissipation terms 
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which are defined in Fourier space as 







fk = 


{ 






. /fc = — otherwise; 


^fe = 








-li,k < kpi, 


7fc = 


1 


0,kpi < k < kp2, 



Here Rj^it) is a uniformly distributed random number in 
the interval (0, 2tt) for each k. We have applied an im- 
plicit difference scheme that keeps the main property of 
this system — conservation of Hamiltonian in the absence 
of pumping and damping. 

The equations were numerically simulated in the 
periodic domain 2tt x 2tt. The size of the grid was 
512 X 512 points. Gravity acceleration g was equal to 
one. Parameters of the damping and pumping in ® 
were the following: kpi = 5, kp2 = 10, kd — 100. Thus 
the inertial interval is equal to a decade. 

In the simulations we paid special attention to the 
problems which could "damage" the calculations. First 
of all, it is the "bottle neck" phenomenon which was stud- 
ied in the paper |2^ . The effect consists in accumulation 
of energy in /c-space ahead of the dissipation region, if the 
dissipation is too large. This effect is very fast, but can 
be effectively suppressed by a proper choice of damping 
value 72 in dissipation © in the case of moderate pump- 
ing values Fq. On the other hand the Fq value should 
not be too small to secure wave cascade on the discrete 
grid. For the case of capillary waves it was examined in 
detail in The second problem is the accumulation 

of " condensate" in low wave numbers due to inverse cas- 
cade. Buildup of condensate can be overcome by simple 
adaptive damping in the small wave numbers. 

After some time the system reaches the stationary 
state, where the balance between pumping and damping 
takes place. In this state an important parameter is the 
ratio of nonlinear energy to the linear one {Hi + H2)/ Ho- 

For example, for the external force Fq = 2x 10~^, 71 = 
1 X 10^'^, 72 — 665 the level of nonlinearity was equal to 
{H1+H2)/ Hq ~ 3 X 10^'^. The Hamiltonian as a function 
of time is shown in Fig. ^ 

The averaged spectrum of surface elevation < \rik\'^ > 
appears to be power-like in the essential part of inertial 
interval, where the influence of pumping and damping 
was small. This spectrum is shown in Fig. [3 

One can estimate the exponent of the spectrum. Com- 
pensated spectra (e.g. multipHed by k^, to get horizontal 
line) are shown in Fig. 13 This figure seems to be the ev- 
idence that the Kolmogorov spectrum predicted by weak 
turbulence theory fits better the results of the numerical 
experiment. 

The quality of the result (closeness to the < \r]k\'^ >~ 
A;~^/^) crucially depends on the width of the inertial in- 



o 

1 

CD 
I 




5000 10000 15000 20000 25000 30000 
Time 



FIG. 1: Hamiltonian as a function of time. Quasi stationary 
state is formed at f ~ 5000. 
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FIG. 2: Averaged spectrum of surface elevation < |77fep >. 
Line ^ kJ^"^ is also shown. 
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FIG. 3: Compensated waves spectra. Solid line — weak tur- 
bulent theory; dashed line — Phillips theory. Only one decade 
in k is shown. 
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terval. In our previous work |28l| similar simulations were 
performed on the grid 256 x 256. Weak turbulent spec- 
trum is clearly seen on the grid 512 x 512, can be hinted 
at the grid 256 x 256, and is almost invisible on the grid 
128 X 128. This difference is demonstrated in Fig. 0] 
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FIG. 4: Compensated spectra for the different grids. 

In the end we would like to mention that for the dif- 
ferent pumping level surface elevation spectra differ only 
due to the different energy flux value P in © , as clearly 
seen in Fig. [S] 
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FIG. 5: Compensated spectra for the different pumping level. 
Solid line — Fo = 3 x 10"", dashed line — Fo = 2 x 10"". 
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